Bell's inequality and universal quantum gates in a cold atom chiral fermionic p-wave superfluid 
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We propose and analyze a probabilistic scheme to entangle two spatially separated topological qubits in a 
Px + ipy superfluid using controlled collisions between atoms in movable dipole traps and unpaired atoms 
inside vortex cores in the superfluid. We discuss how to test the violation of Bell's inequality with the generated 
entanglement. A set of universal quantum gates is shown to be implementable deterministically using the 
entanglement despite the fact that the entangled states can only be created probabilistically. 
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Introduction: Topological quantum computation affords 
the amazing possibility that qubits and quantum gates may 
be realized using only the topological degrees of freedom of 
a system lit]. Since these degrees of freedom, by definition, 
are insensitive to weak local perturbations, the resulting com- 
putational architecture should be free of environmental deco- 
herence, a major stumbling block to quantum computation. 
In a class of topological systems, the requisite (non-Abelian) 
statistical properties 12, IH are provided by the presence of 
Majorana fermion excitations described by the self-hermitian 
operators 7^ = 7. These excitations have been shown to oc- 
cur naturally at the cores of vortices in a 2D spinless + ipy 
superfluid or superconductor where the interacting 

fermions are described by the many body Pfaffian wavefunc- 
tion |i2J. (It seems likely, but remains to be verified, that this 
wavefunction also describes the essential physics of the fill- 
ing fraction v = b/2 fractional quantum Hall (FQH) system 

^). It is encouraging that the spinless px + ipy superfluid 
of fermionic cold atoms is potentially realizable in an opti- 
cal trap tuned close to a p-wave Feshbach resonance iSSMl- 
Our current work establishes the possibility of testing Bell's 
inequality in a cold atom p-wave fermionic superfluid on the 
way to eventual universal topological quantum computation 
using vortices in such a system. 

In a + ipy superfluid, one can define a topological qubit 
using a group of four vortices. Since the states of the qubit are 
associated with the composite states of the four spatially sepa- 
rated Majorana fermion excitations, they are immune to local 
environmental errors. One can implement some single-qubit 
gates by adiabatically moving (braiding) one vortex around 
another within the same vortex complex defining the qubit. 
Since the associated unitary transformations are purely statis- 
tical, there is, in principle, no error incurred in these gating 
operations. However, it is well known [IQJ that such a braid 
operation of one vortex from one qubit around another from a 
different qubit fails to provide a two-qubit gate: the topolog- 
ical braiding operations allowed in a p^ + ipy superfluid, as 
in its FQH Pfaffian counterpart, are not computationally suffi- 
cient. 

The principal reason why a p^; + ipy superfluid is not com- 
putationally universal is that two qubits cannot be entangled 
using only the topological braiding operations. Any compos- 
ite state of the two qubits, accessible by braiding one excita- 



tion around another, can always be written as a product of the 
states of the individual qubits [T^. Therefore, in light of its 
experimental relevance, it is important to examine the problem 
of creating quantum entanglement in a p^; + ipy superfluid via 
some other, possibly non-topological, means (without incur- 
ring too much error) which, coupled with the available braid- 
ing transformations, may lead to universal quantum computa- 
tion. This is all the more important because the other, more 
exotic, non-Abelian topological states, e.g. the SU(2) Read- 
Rezayi state ifTTIl . which can support universal computation 
via only the topologically protected operations ifTil. Fisl] . are 
presently much beyond experimental reach. In the 5/2 FQH 
state, non-topological interference of charge-carryi ng q uasi- 
particle currents along different trajectories i fiol Il4l I1S I was 
proposed to entangle qubits. Such an approach is not suitable 
for the superfluid, because the non-Abelian excitations here 
are vortices, which do not carry electric charge. 

In this Letter, we show how to entangle two spatially sepa- 
rated topological vortex qubits in a cold atom px + ipy super- 
fluid by using two other, movable, external dipole traps. (The 
two-state systems formed by the atoms in the movable, exter- 
nal traps will be referred as ''flying qubit"). Controlled cold 
collisions between an atom in the dipole trap and an atom at 
the vortex qubit yield entanglement between the flying qubit 
and the vortex qubit. Subsequently, a measurement on a sys- 
tem comprising two flying qubits, entangled with two differ- 
ent vortex qubits, collapses the two vortex qubits on an entan- 
gled state. We show how to test the violation of Bell's inequal- 
ity with the obtained entanglement. Finally, we show how to 
deterministically implement a set of universal quantum gates 
using the entangled state, although the entanglement among 
the vortex qubits itself can only be generated with a 50% suc- 
cess probability. It is important to mention that the entangle- 
ment can be generated and purified off-line, and so the non- 
topological nature of the corresponding operations does not 
degrade the topological quantum computation. 

Topological qubit and flying qubit: Consider a quasi-two 
dimensional (xy plane) Px + ipy superfluid of spin-polarized 
atoms where vortices in the superfluid can be gen- 

erated through rotation or external laser fields. For each vor- 
tex, there exists a zero energy state that supports a Majorana 
fermion mode 7 ilia. Two Majorana fermion states in 
two vortices can be combined to create an ordinary fermionic 
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Figure 1: (a) Single qubit flip gate R = —ia^- (b) Single qubit phase 
gate A (7r/2) — diag(l, i). (c) Hadamard gate H = ( ^ \ / ' 



state c = (71 + 172) /2. Therefore, a natural definition of a 
vortex qubit may be given in terms of the unoccupied, |0), 
or occupied, |1) — c' |0), states of two Majorana vortices. 
Here the occupied state |1) contains an unpaired Fermi atom 
inside the cores of the vortex pair ("unpaired" as opposed to 
"paired" in a cooper pair in the superfluid). However, such a 
definition does not allow the superposition of the basis states, 
i.e., the states, (|0) ± |l))/\/2, do not exist because they vio- 
late the conservation of the total topological charge (the super- 
fluid condensate conserves the fermion number modulo 2). To 
overcome this difficulty, a topological vortex qubit is defined 
through two pairs of vortices, i.e., with the states |0)y = |00) 
(the two vortex pairs, (1,2) and (3,4), are both unoccupied), 
and \l)y = |11) (the two vortex pairs are both occupied). 
The superposition states, (|0)y ± V2, are now allowed. 

Note also that these two states do not mix, via any unitary 
braiding operations, with the other two states of the four- 
vortex complex, 1 10), 1 01). Various intra- and inter-pair vortex 
braiding operations within a single qubit give rise to various 
single-qubit gates lfl4l[l6ll (e.g. qubit-flip gate R, phase gate 
A(7r/2) and the Hadamard gate H) as depicted schematically 
in Fig. [T] Finally, the state of the vortex qubit can be read out 
in the {|0)y , |l)y} basis by fusing the vortices pairwise and 
detecting the number of unpaired atoms in the core [6]. 

The flying qubit is constructed using an atom trapped in the 
ground state of a movable optical dipole trap which is itself 
formed by overlapping two identical laser beam traps. One 
laser beam trap can then be adiabatically moved out to split the 
composite trap into two traps, L, R, see Fig.|2a). This yields 
a superposition state for the atom, (|01)^^ + |10)^^) /\/2- 
Here L and R denote the left and the right traps, respectively. 
Now, concentrating on the left trap only, one can define the 



two states of the flying qubit, \0) p — |01) 
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Note that the two states of the qubit are distinguished by the 
absence (|0)^) and the presence (|l)p) of the atom in the left 
dipole trap. 

Entanglement between two topological qubits: As is well 
known [IQil , two topological qubits cannot be entangled by 
braiding one vortex from one qubit around another from the 
second qubit. However, using the flying qubits as auxiliary 
degrees of freedom, one can generate entangled states be- 
tween the two qubits. The basic idea of the entanglement 
generation is illustrated in Fig. |2b, c). Initially, a vortex 
qubit, V, is prepared in the state |0)y. A Hadamard gate 
H is applied to the qubit that transfers the state to 1 0) ^ — 
(|0)y + |l)v) /\/2- By splitting a composite dipole trap in 
two parts (Fig. |2|a)), the flying qubit F is prepared in the 



state \ip)p = {\fS)p + \l)p) 1^/2. The flying qubit is then 
moved near to the vortex qubit (Fig. |2|b)) so that the trapped 
atom can collide with the unpaired fermi atom, if any, in the 
vortex pair As shown below, such a collision process yields a 
controlled phase gate, CP(7r) = exp {iirnvnF), between the 
flying qubit and the vortex qubit, where ny = 0, 1 is the num- 
ber of unpaired fermi atom in the vortex pair and np = 0, 1 
is the number of atom in the flying qubit. It is easy to see that 
the gate CP(7r) gives rise to the transformation. 
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(|o)v + |i)v) + |i>i.(|o)y-|i)^.)]A 



which can be transferred to an entangled state 
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= (|0)^|0), 



|l>Fll>y 



/V2 



(1) 



between the flying qubit and the vortex qubit by applying a 
Hadamard gate on the vortex qubit. 

Two vortex qubits can be entangled by a projection mea- 
surement on the flying qubits of two entangled states 1$)^^^^ 
and \^) p^y^. The dipole traps of the two flying qubits are 
spatially merged and the atom number is measured through 
fluorescence signals (Fig.|2c)). From the combined state, 

\^)f,vA'^)f,v, = ^(|00)^,^J00)^^^^ + |11)^^^J11)^^^^ 
+ |01)f^,^j01),.^^^ + |10)^^^J10) 

where lOO)^^^^ ~ \^) Fi \^) F2 ^^^^ '•° deduce the 

probabilities for the three possible outcomes: one atom (50%), 
zero atom (25%), two atoms (25%). In the last two cases, 
the states of the vortex qubits are projected to \^)yi |0)^ and 
IDy, ID 



respectively, and are not entangled. Therefore, 
in these cases the above procedure for creating the entangled 
states, |*I')^jyj^ and |$)p^y^, need to be repeated. However, in 
the case where the measurement produces one atom, the quan- 
tum state of the two qubits is projected to the entangled state 

(|0)yi |l)v2 \^)vi \^)v2) /"V^' which can be transferred to 
the expected entangled state 



= (|0>yi |0)^ 



+ |1)^,J1)^,J/V2 (2) 



using simple qubit-flip gates. Note that the above entangled 
state can only be created with a 50% success probability. For 
later use, the gate representing the generation of entanglement 
is denoted as "EG". 
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Figure 2: (a) Construction of the flying qubit by sphtting a composite 
dipole trap into two traps, (b) Reahzation of the gate CP(7r) by con- 
trolled collisions of atoms, (c) Two flying qubits are merged into one 
and the number of atoms is measured through fluorescence signals to 
create entanglement between two topological qubits (see text). 
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The remaining problem for the entanglement generation is 
how to realize the controlled phase gate, CP(7r), between the 
flying qubit and the vortex qubit. In Fig|2]3, the center of the 
dipole trap, Fq (t) = Zq (t) (with the core of vortex 1 as 
origin) is adiabatically brought from a distance doe^ above 
the z = plane, where the wavepackets of the dipole trapped 
atom and the unpaired atom in the vortex pair (1,2) do not 
overlap, to a distance zero, where they do. The collision 
phases between the atoms are dynamic phases. They are dif- 
ferent for different quantum states of flying and vortex qubits 
with different total energy. 



E = Ef (i) + Ev {]) + AE, (z, j) , 



(3) 



where, i,j ~ 0,1 correspond to the quantum states 
|0) and respectively. Ep (0) = and Ep (l) = 
Jd^ja* (r ~ fo {t))[~h^Vy2mp + Vpir~ fp (t))] 
a(f — fo {t)) + Eg are the energies of the flying qubit in 
the states |0) and respectively. Vp — i^o (t)) is the 
harmonic potential of the dipole trap, and a (f — Fo {t)) is the 
ground state wavefunction of the dipole trapped atom with 
mass mp- Eg is the interaction energy between the dipole 
trapped atom and the paired BCS condensate. The second 
term Ey (j) corresponds to the energy of the fermionic 
state in the vortex cores near the dipole trap. Because these 
states are the solutions of the Bogoliubov-de Gennes (BdG) 
equations with eigenvalue zero, Ey (j) — for j = 0, 1 HQ]. 

The last term describes the collision energy UtIi between 
dipole trapped atoms and unpaired fermi atoms, and is non- 
zero only if both the flying qubit and the vortex qubit are in 
the occupied state, 

Ai?e(l,l) = | j d'v\a{v-Vo{t) fSnv{v). (4) 

Here g is the collision interaction strength and Sny (?) de- 
notes the changes of the atom density from the BCS conden- 
sate density. Using the standard harmonic trap wavefunction 
for a (r — fo {t)) and the wavefunction for the zero-energy 
mode obtained from the solution of the BdG equations, we 
find 



ASc (1, 1) = hn± cxp {-zl (t) /a^) 



(5) 



with a D and ay the oscillation lengths 



where = ajj-{^ 
for harmonic trapping potentials along the z direction of the 
dipole trap and the superfluid, respectively. h^l± is the char- 
acteristic energy scale for the collision interaction, which is 
determined by the collision interaction strength g as well as 
the overlap between the wavefunctions of the dipole trapped 
atoms and unpaired fermions in the vortex pair (1,2). Note 
that the occupied zero energy fermionic state of the vortex 
pair (1,2) is as likely to contain extra atoms (corresponds to 
il+) as to miss atoms (corresponds to f2_) compared to the 
BCS condensate, which yields $7+ = — ri_. 

The state-dependent energy Q yields a state-dependent dy- 
namic phase 



where (fp (i) = Jl^Ep {i)dt is the single qubit phase 
that can be incorporated in the definition of the flying 
qubit, the collision phase (/)c(i,j) = J__^ AE^ii, j) dt, 
and =FT denote the time when the dipole trap center fo (t) 
moves from and back to the initial place doe^- As- 
suming that Tq {t) varies adiabatically as zo{t)/dQ = 
T] (exp (i^/r^) — l) / (l + 7/ cxp (t^/r^)) with the parame- 
ter rj ~ exp (— rf /t^), the controlled collision phase can be 
written as. 



{l,l) — n±TrJ exp 



-Try- 



dt (7) 



where T — d^/a? and time in the above integration has been 
scaled by t^. With a set of parameters for ^Li, ap, ~ ay — 
OAfj,m, dg = lOap) = 4/im, = = 3.57/17, r = lOr^, 
s-wave scattering length ^ 53nm, and the vortex core size 
^ ^ l/im, we estimate fl± ^ ±2tt x 6.6kHz, r ~ 0.86ms 
and the phase 9 = zLtt (i.e., exp {iff) = —1). Therefore the 
controlled phase gate CP(7r) can be realized. 

Violation of Bell's inequality: The entangled state l^)^^^-^ 
between two remote vortex qubits can be used to test the vi- 
olation of the CHSH inequality, a variant of the Bell's in- 
equality 1 18]. Violation of the CHSH inequality would estab- 
lish the quantum non-locality between the two vortex qubits. 
A schematic diagram of this test is given in Fig. [3] The 
test requires to measure the vortex qubits along four differ- 
ent directions: Ai = aY^ (g) A2 = a^^ « Bi ^ 



-/^i ® {aY^ + a^^) /V2, B2 = A ® {aY^ - a^^) /^/2. 
After the measurements, two parties at V\ and V2 need to 
communicate their results through classical channel. After re- 
peated measurements, the statistical average L = {AiBi) + 
{A2B2) + {A2B1) — {A1B2) is evaluated. The quantum non- 
locality of the entangled state yields L = 2v^, which violates 
the CHSH inequality for local realism, L < 2 

It is easy to convince oneself that the above four mea- 
surements correspond to measuring the two vortex qubits 
in four different bases which are eigenstates of the re- 
spective operators, Ai : Vi on {|0)y^ , |l)y^}; A2:Vi 



on 



{(|0> 



Vi 



1 



Vi 



/V2,{\0)y^-\l)y^)/V2}; B, : V2 
on {a\0)y^+b\l)y^,b\0)y^-a\l)yj- B2 : V2 on 
{a|0)^,^ -6|1)^^,6|0)^^ +a|l)yj, where a = cos (^8), 
b ~ sin(7r/8). In the experiment, Ai is a fusion mea- 
surement of the number of unpaired atoms in the vortices 
[6]. Measurements A2, Bi, and B2 can be implemented by 
first applying suitable single-qubit operations to the qubits 



(a) 
(b) 



A; 



EG 



Classical Communication Channel 



- Sj - I -\'-{;r/2)^W]—A{7r/4)-\l{]- A, 



(6) 



Figure 3: (a) Testing the violation of the CHSH inequality, (b) The 
realization of the Bi measurement in (a). 
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to transfer their measurement bases to {|0)y , follow- 
ing by fusion measurement Ai. The corresponding single- 
qubit operations are A2: H; Bi. HA (e*''/'*) HA'^ {tt/2); B2: 
HA {e"/^) HA^ (-vr/2), where A (e''^/'') = diag(l, e"/"*) 
is a single qubit phase gate. A (e*'^/'') cannot be implemented 
through topologically protected braiding operations and its re- 
alization is discussed in the next section. 

Universal quantum gates: It is well known that a set of 
quantum gates [10, 14] 

H, A(e*-/4^, A (a,) (8) 

are sufficient to simulate any quantum circuit, where A (cr^) — 
diag(l, 1, 1, — 1) is the two-qubit controlled phase gate be- 
tween two vortex qubits. Among these three gates, the 
Hadamard gate H can be implemented using the topological 
braiding operations. The single-qubit phase gate A (e"/*) 
can be realized by bringing two vortices 1 and 2 of a 
qubit close together il4tl . The tunneling between two vor- 
tices induces an energy splitting /S.Ez [t) between the state 
|0)y and which yields a relative dynamic phase exp 

{i A_E^ {t) dtj between these two states, with 2Tp as 
the total tunneling period. A properly chosen tunneling pro- 
cess with /^^ AEz (t) dt = 7r/4 yields the gate A (e*^/"*). 

A controlled phase gate A (az) between two arbitrary vor- 
tex qubits can be realized deterministically provided one has 
been able to create the entangled state I'i) between two vor- 
tex qubits. Considering two vortex qubits G and Q (with 
the constituent vortices Gi, G2, Qi, Q2 etc.), we note that 
Agq (ct^) = Ag {tt/2)Aq (7r/2)exp(i7r7Gi7G2 7Qi7Q2/4), 
where the last term involves interaction among four vortices. 
The requirement of a four-vortex interaction is indeed the 
reason why the two-qubit gate cannot be implemented using 
braiding operations which can lead to only two-vortex (sta- 
tistical) interactions. The four- vortex operator can be imple- 
mented using one additional vortex pair (7W1 , 711^2 ) (initially 
prepared in state |0)) by noting that 1 19], 

exp (z^7Gi7G2 7Qi7Q2/4) = 2U^,P^^^ p!,^\ (9) 

where = (It ilQ,lw,)/2 and Pi^^ = 

{I ± 7Gi7G2 7Q2 71^1)72 are non-destructive measure- 
ments which project the state of the vortices to the eigen- 
states of the operators —ijQiJWi and 7Gi 7G2 7(32 7^1 ■ 
Ufj,„ are corresponding braiding operations for dif- 
ferent measurement results {/ii^}, — = 
e'^7Qi7»'2/4, = iAG(j)AQ(i)e'''^«i'^»'2/4^ 
U-+ = iAG(«) Aq (i)e~''''«i'''^2/4. Here e''""^^""^^/'^ 

is just the exchange of the vortices 7q ^ and 7VK2 ■ 
(2) 

Pj_ ' can be realized via a basis transformation method. 
We exchange the vortices 7Qj and -fwi to transfer 

two eigenstates of -ijQ,-/Wi to ||00)g^y, |11)qvv'} or 

||10)q^, |01)q)^|, depending on the total topological 
charge of the four vortices jq-^, -fQ^, jwi and 71^2- We then 



apply a fusion measurement on the vortex pair {jwi , IW2 ) 
to determine whether the state is |0)^y or which corre- 

spond to the eigenvalues +1 or —1 of the projection measure- 
ments Pj. . After the fusion measurement, the vortex pair 
(7Wi> IW2) is recreated in the state |0)^;;^. If the result of 
the fusion measurement is the state this state is recov- 

ered by applying a single-qubit flip operator R. Vortices 7Qj 
and are exchanged again to transfer the states back to 
the eigenstates of ^i^Q^^Wi- With this basis transformation 
method, the projection measurement Pj. can be performed 
non-destructively. 

However, such basis transformation method does not work 
for the measurements P^*-* because they involve eigenstate 
measurement of four vortices. Recent work |10] showed 
mathematically that p'^^ can be realized deterministically us- 
ing the auxiliary entangled state \^), for which we provide 
a prescription in this Letter, coupled with the braiding opera- 
tions and the fusion measurements. Here we refer the math- 
ematical details of this measurement to Ref. ifioll . Note that 
the measurement p'^^ can be deterministically implemented, 
although 1^1') in our scheme can only be generated with a 50% 
success probability. This is because j^I') can be prepared using 
off-line procedures that are not involved in the measurement 
process. In addition, pairs with non-perfect entanglement can 
be purified to pairs of nearly perfect entanglement through off- 
line purification processes. Therefore, the controlled phase 
gate A [az] can be implemented with a high accuracy because 
the remaining processes only involve the braiding operations 
and the fusion measurements. 

In summary, we proposed and analyzed a scheme to gen- 
erate entanglement between two topological vortex qubits in 
a px + ipy atomic superfluid with the assistance of external 
flying qubits. The entanglement can be created and purified 
off-line and therefore, in spite of being a non-topological pro- 
cess, does not degrade the actual quantum computation which 
continues to use the topologically protected braiding opera- 
tions. We showed how to test the violation of Bell's inequality 
using the obtained entanglement. Finally, we showed how to 
deterministically implement a set of universal quantum gates 
in the chiral p-wave superfluid, which had hitherto remained 
a major conceptual problem, using the entanglement created 
between two topological qubits. 

We thank A. Kitaev for helpful discussion. This work is 
supported by LPS-NSA and ARO-DTO. 
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